
% Parametros
intervalos = 5;
alfa = 0.05;

% Genero el vector de numeros pseudo-aleatorios.
load 'numeros.txt', numeros;
cant_numeros = length(numeros);

% Genero las marcas de clase 'x' en funcion de la cantidad de intervalos.
vmin = min(numeros);
vmax = max(numeros);
ancho = (vmax - vmin) / intervalos;
x = vmin + ancho/2 : ancho : vmax - ancho/2;

% Los clasifico en intervalos y obtengo la frecuencia absoluta 'f'
[f, x] = hist(numeros, x);

% Valor teorico esperado de las frecuencias de cada intervalo.
teorico = cant_numeros / intervalos
f

% Calculo el estadistico chi
chi = 0;
for i=1:intervalos
	chi = chi + (f(i) - teorico )^2;
end
chi = chi / teorico

% Obtengo el valor critico por tablas.
chi_tabla = chisquare_inv(1 - alfa, intervalos - 1)


